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Effects of H 2 coating of grains on depletion of molecular species 
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ABSTRACT 

Physical conditions in dense and cold regions of interstellar clouds favour the formation of ice 
mantles on the surfaces of interstellar grains. It is predicted that most of the gaseous species 
heavier than H2 or He will adsorb onto the grains and will disappear from the gas-phase, 
changing its chemistry, within ~ 10 9 /nn years. Nonetheless, many molecules in molecular 
clouds are not completely depleted in timescales of 10 5 yr. Several speculative mechanisms 
have been proposed to explain why molecules stay in the gas phase, but up to now none are 
fully convincing. At the same time, these mechanisms are not mutually exclusive and we 
can still explore the effects of other possible processes. We speculate on the consequences 
of H2 coating of grains on the evaporation rates of adsorbed species. More experiments and 
simulations are needed to calculate the evaporation rate i? ovap (X-H2). 
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1 INTRODUCTION 

Icy mantles on dust grains have been l ong detected in cold a nd den- 
se regions of the interstellar medium dWilliams et al.|[l992t) . These 
mantles are the result of the accumulation of heavy molecules 
on grain surfaces due to the low temperature and high density 
conditions of the environment. The collisions between the dust 
grains and the neutral component of the gas lead to the reten- 
tion of the gas phase particles on the grain surface, because 
their thermal energy is lower than the typical adsorption energies 
and because th e excess kinetic energy can be transferred rapidly 
to the surface l|Herbst & Millar! 1 19911) . Calculations have shown 
dLeitch-Devlin & Williams! 1 9851) that ~ 30-100% of the collisions 
result in retention, depending on the species and the nature of the 
surface. Moreover, the widespread detection of molecular mantles 
in dense regions of molecular clouds indicates a high efficiency of 
the sticking of heavy atoms and molecules on the surfaces of cold 
dust grains. 

The main observed component of the mantles is water ice 
dTielens & Hagenll 19821 : Iwhittet et al.ll 19981) , but it is thought that 
most of the gaseous species heavier than H2 or He will adsorb 
onto the grains and will disappear from the gas-phase, changing 
its chemistry, in a relative ly short "freeze-out" time. The fre eze-out 
rate is ~ 10 9 /riH years dLeitch-Devlin & Williamsll 19851) . where 
tih = n(H) + 2n(H2) is the total hydrogen nucleon number den- 
sity. For densities as low as ~ 10 4 cm~ 3 , the freeze-out timescale 
is much less than the expected lifetime of a typical molecular cloud 
f~ 10 7 yr), so it would be expected that, contrary to the obser- 
vational result, the majority of ob servations wou ld not show evi- 
dence for heavy gas-phase species. Ilglesiasl ( Tl977l) showed that for 
non-zero sticking efficiency, when no desorption processes act to 



remove mantle material, the accretion of gas-phase species on to 
the dust dominates the chemical evolution and the total freeze-out 
of gas is inevitable. 

Thus, some mechanism is preventing the growth of mantles in 
grains and maintains an equilibrium between molecules in the gas 
and solid phases by producing the desorption of molecular species 
from the surface of dust grains. The finding of a convincing des- 
orption mechanism that explains the observed molecular gas -phase 
abundances has proved to be rather elusive in the last 20 years. It 
is not yet precisely known how this desorption takes place and sev- 
eral different continuous desorption processe s have been propose d, 
either working simultaneously or separately dRoberts et alj|2007h : 

- Classical thermal evaporation dLegeretal.lll985l) . which is 
negligible at the temperatures of dense clouds, ~ 10 K 

- desorption by direct impact of cosmic rays (or X-rays) onto 
grains, which causes local heating that lead to the evaporation of 
weakly bound molecules, such as CO dWillacv & Williams|[l993l : 



iLeger et al.lll985l ; lHasegawa & Herbstll 19931) 



cosmic ray-induced photo-desoiption dHartquist & Williams! 
Il99d : lwillacv etal.ll 1994) 

- exothermic mantle reactions, such as the formation of H2, 
which can result in a local "hotspot" heating of the grain mantle 
dWillacv et al.lll994l : lDulev & Williams! 19931) 

- g rain-grain collisions or chemical explosions dLanger et al.l 
120001); grain-grain collision or sputtering in interstellar shocks 
dTielens et anil 9941) 



- if the cloud is composed of clumpy molecular gas layers, in- 
terstellar UV radiation can penetrate two or three times deeper than 
in the case of a homogeneous cloud and can produce orders of mag- 
nitude more photo-io nisation and photo-desorption eve n at modest 
extinctions, Ay ~ 2 dBoisse1ll99ol : lBetheil et ai]l2007t) . 
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The feedback from the star formation processes, which 
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eventually disrupt the dense cores from which stars form, will 
also limit the growth of ice mantles or force the return of 
the adsorbed species to the gas phase, in a sporadic or inter- 
mittent manner, through a variety of mechanisms: by sublima- 
tion of icy mantles from warm grains near brig ht stars in a 
more or less gradual way ( Viti & Williamsl Il999t ISchoier et al.l 
120021 : iRodgers & Charnlevl 120031 : iGarrod et al.l 120081) . by de- 
structio n or sputtering of gr a in cores and mantles by passing 
shocks dCharnlev et alj Il988l: iFlower & Pineau des Foretsl Il994l : 



iBergin et alJI 19981 : Ivan Dishoeck & Blakdll998l ), or by ice subli- 
mation due to UV r adiation originating in shocks dViti et af]|2003l : 
IChristie et all201lh . 

All the mechanisms above are not mutually exclusive, but no 
proposed non-thermal mechanism convincingly explains the pre- 
sence of molecules such as CO in dark clouds or it is widely 
accepted as a reliable desorption mechanism (e.g. , iFrerking et al.1 
ll982l : lDuvert et al.lll986l : iDickman & Herbstlll990l) . In this paper, 
we propose to explore another possible mechanism: the conse- 
quences of the variation of the adsorption energies due to the dif- 
ferent properties of the substrate and, in particular, the effect that 
the coating of dust grains by molecular hydrogen has on the evapo- 
ration rates of adsorbed molecules. 



considerably reduced when compared with the ones for silicate or 
H2O ice surfaces. For instance, the adsorption energy for H2 on 
H 2 ice, E D = 350 K, while on H 2 ice, E D = 100 K. The adsorp- 
tion energies of other accreting heavy molecules behave in a similar 
way. Once En's are reduced to 1/5 of the silicate surface values, 
evaporation timescales would be shorter than freezing timescales, 



2.1 Coating of grains by H 2 

By definition, H2 is the most abundant gas-phase species in mole- 
cular clouds and several authors have long susp ected that it could 
be an important compone nt of the grain surfaces dfielens & Hagenl 
ll982tlGoversetalJll98Cl) . At the same time, the adsorption ener- 
gies of H2 indicate that it is a rather volatile molecule. Table [TJ 
shows that the evaporation time for H 2 on H2O ice can go from a 
few seconds at a grain surface temperature of 15 K to 0.5 years for 
10 K. 

The surface coverage of H2 molecules will depend on the ba- 
lance between their evaporation and accretion rates. The number 
of accreting H2 molecu les per unit of time can be calculated from 
iTielens &Hagenld 19821) 



2 EFFECTS OF H 2 COATING OF GRAINS ON 
ADSORPTION ENERGIES 

Once a molecular species sticks to a grain, the time it stays on 
the surface before it evaporates depends on the surface tempera- 
ture and the binding (adsorbing) energy of the adsorbate species 
on a given subs t rate. T he desorption rate, in s _1 , is described by 
lHasegawa et alJdl992l ) 



-Rcvap = vq exp(-E D /T) 



(1) 



where T is the surface temperature of the grain, Ed the adsorption 
energy in K, and vq is the characteristic vibrating frequency for the 
adsorbed species, uo = (2n s Eo /7r 2 m) 1,/2 , where n a is the surface 
density of sites (« 10 15 cm -2 ) and m is the mass of the adsorbing 
particle. Typically vq ~ 10 12 -10 13 s _1 . The adsorption energies, 
Ed, of molecules onto the grains depends on the properties of the 
substrate. Table Q] shows the adsoiption energy, in K, of several 
fundamental species on different substrates that can be found in 
molecular clouds. TableQ]also shows the evaporation time expected 
for each adsorption energy for two relatively close grain surface 
temperatures, Td us t = 10 and 15 K. It is easily visible the large 
differences between the different molecular species and substrates 
and the influence of temperature. 

When modelling the accretion of gas-phase molecular 
species onto the grain surfaces, many adsorption energies are 
needed to describe the complex interacting system thus created 
dCuppen & H erbst 2007]). Ideally, we would need to know the ad- 
sorption energies of the species onto the carbonaceous substrate, 
but also the energies between the species and ice (or different ice 
structures) and with other adsorbed species. U nfortunately, many 
of the se energies are poorly known, if at all (see lCuppen & HerbsJ 
120071 and references therein). The lower half of Table [T] also 
shows some estimates of the adsorption energies of several molec- 
ular species when the substrate is H2. Most of these values were 
calculated assuming, for demonstration purposes and following 
lAllen & Robinsonl dl977l) . a scaling value of 1/5 of the adopted 
standard value of Ed - The adsorption energies for H2 surfaces are 



Rac = s d {VH 2 )n(R2)S = 0.10193 



n(H 2 



10 4 c 



(2) 



where n(H2) is the abundance of gas-phase H2 molecules, (Vk 2 ) 
is the mean velocity of H2 molecules (~ 1.026 x 10 4 v / T\ cm s _1 , 
where Tk is the gas kinetic temperature), Sd is the grain geomet- 
rical cross-section (~ 3.14 x 10~ 10 cm -2 ), and S is the sticking 
coefficient, which we will assume equal to 1. In practical terms, for 
a dense gas density, ~ 10 4 cm -3 , a H2 molecule should hit a grain 
every 10 s. 

The number of H2 molecules that evaporate from the surface 
of a grain per unit of time can be expressed as 



Reva P = iV(H 2 ) v exp(-£ D (H 2 )/T) 



(3) 



where Af(H 2 ) is the number of H2 molecules on the grain surface. 

If the number of grain sites, N s = 10 6 , the accretion time is 
~ 10 7 s or ~ 0.3 yr. For a dust temperature of 10 K, and _Ed(H 2 — 
H 2 0) = 450 K, the evaporation time is ~ 1.47 x 10 7 sor~ 0.5 
yr. The time needed for the system to reach steady state can be 
calculated as 



riiV(H 2 ) 
dt 



Rac Re 



and the number of H2 molecules on the grain surface will be 



JV(H 2 ) = 1.5 x 10 6 



n(H 2 



10 4 cm- 2 



(4) 



(5) 



For a number of grain sites, N s = 1-3x10 , and n(H 2 ) > 
10 4 cm~ 3 , all the exposed sites will be occupied by H 2 molecules 
in ~ 0.5 yr at nn > 10 4 cm -3 dTielens & Hagenll 19821) . The evap- 
oration of H2 from an H2 ice surface is very fast, leaving no possi- 
bility of multiple layers of H2 ice. 

Given the possibility that grains in molecular clouds are co- 
vered by a layer of H2 molecules and, as Table Q] shows, this can 
greatly reduce the adsorption energies of heavy molecules, it is re- 
levant to explore how much do the changes in the adsorption ener- 
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Table 1. Adsorption energies and evaporation timescales for different grain surfaces. 



c : TT 

Species H 


H2 


H2 


He 


O 


CO 


CO 


CO 


Surface H2O ice 


H2 O ice 


H2 ice 


Silicate 


H2O ice 


graphite 


H2 O ice 


CO ice 


E D (K) 350 a 


450 a ' 6 


100 c 


100 a 


800 a 


1208 d 


1740 c 


960 d 


Tkust = 10 K 
















^cvap 530 S 


0.46 yr 


2 x 10~ 8 s 


3 x 10~ 8 s 


2 x 10 15 yr 


9 x 10 32 yr 


9 x 10 55 yr 


2 x 10 22 yr 


Td us t = 15 K 
















^cvap 5 X 10 S 


4.5 s 


7 x lO" 10 s 


1 x 10" 9 s 


4 x 10 3 yr 


3 x 10 15 yr 


6 x 10 30 yr 


2 x 10 8 yr 


On H2 coated grains 
















Species CN 


NH 3 


He 


O 


N 2 


CO 


CO2 




Surface H2 ice 


H2 ice 


Silicate 


H2 ice 


H2 ice 


H 2 ice 


H2 ice 




E D e (K) 302 


615 


100 


160 


340 


348 


572 




Tkust = 10 K 
















^cvap 24 S 


1.7 x 10 7 yr 


3 x 10~ 8 s 


2 x 10~ 6 s 


18 min 


40 min 


3.8 x 10 5 yr 





a from lTielens & Allamandolal <1987l) . 

6 450 K was calculated bv lTielens & Allamandol3 ll987l) based on lHollenbach & Salpetej jl970h . 
c from lSandford & Allamandolal l ll993al) . 
d from lAUen & Robinsonl f 19771) . 

e Values estimated using a scaling factor of 1/5 of the adopte d standard value of the a dsorption energy, E D jAllen & Robinsoijl 19771). The 
origina l standard values of E p were: CN on graphite 1510 K fallen & Robinsonll 19771); NH 3 on NH 3 ice, 3075 K ISandror^&Allamandor3 
Il993ah: O on H 2 Q ice, 800 K iTielens & Allamandolal 19871) : N 2 on H 2 Q ice, 1700 K iTielens & Allamandolal 19871) : CO on H 2 ice, 1740 K 
ISandford & Allamandolal Il993ah : and COo on H2O ice, 2860 K {Sandford & Alla mandolj|l993 j) . 



gies used in the chemical models change the values of molecular 
depletion in the gas. 



3 RESULTS OF THE CHEMICAL MODELLING 

We used our gas-grain chem i cal m odel based on the gas-grain 
chemistry of lHasegawa et all i 19921) to explore the evolution of 
gas- and solid-phase abundances of species after assuming that the 
grains are coated with H2 molecules. The chemical model con- 
tained 3768 gas and surface reactions involving 503 species, but 
we only took into account adsorption and evaporation of species 
from the grain surfaces, with no other chemical reaction on them 
apart from the formation of H2. We ran several models at differ- 
ent gas densities (2 x 10 4 , 2 x 10 5 , and 2 x 10 6 crrT 3 , at a 
gas and dust temperature of 10 K, and cosmic ray ionisation rate 
£ = 1.3 x 10 -17 s - 1 . All models used initi al elemental 'low metal' 
atomic abundances dRuffle & Herbst 2 000h . 

We assumed a crude approach to simulate the effects of the 
coating of grains by H2 molecules and the subsequent modifica- 
tion of the adsorption energies, Ed- We ran three different sets of 
models: a fiducial one with no modification of Ed, and two mo- 
dels where all Ed, except for H, H 2 and He, were reduced by ei- 
ther a factor of 1/2 or 1/3 , a conservative value stil l far away from 
the scaling factor used bv lAllen & Robinsonl l l 19771) to estimate the 
change in the adsorption energy of molecules landing on H2 -coated 
grains. Figure Q] shows the results for two different number densi- 
ties, 2 x 10 4 and 2 x 10 6 cm , and three different sets of adsorption 
energies: standard values, 1/2 of the standard Ed values, and 1/3 
of the standard Ed values. 

The upper row of Fig.[T]shows the results of the fiducial case. 
For a gas density of 2 x 10 4 cm -3 , the gas-phase molecular abun- 
dances begin to suffer from depletion at a few times 10 yr, at which 
point the gas and grain-phase abundances are approximately simi- 



lar, and most of the molecular species disappear from the gas-phase 
at a few times 10 6 yr. This effect is accelerated by about a factor of 
10 for a higher density, 2 x 10 6 cm -3 . 

The middle row of Fig. [T] shows the results of the case where 
the adsorption energies are half of the standard values. There is still 
some degree of depletion (larger than 90%) for some molecules, at 
a density of 2 x 10 4 cm -3 , but total depletion is never reached. 
Some molecules, such as NH3, do not have so high a depletion and 
the gas-phase abundance remains higher than the solid-phase abun- 
dance at all times. The level of depletion is much larger for higher 
densities, but even in this case there is not an almost total disap- 
pearance of the gas-phase molecules as in the standard case. On 
the other hand, the times when the solid-state abundances become 
larger than the gas-phase abundances do not seem to change that 
much. 

The case where the molecular adsorption energies have been 
reduced by a factor of 3 (bottom row of Fig. [T} shows the more dra- 
matic changes. Many molecules, such as CO, CS (prime candidates 
for depletion) or NH3, show little or insignificant depletion, even at 
high densities, measured as the fraction of molecules in the gas and 
solid phases. This is not the case for CN, which still shows high 
depletion in this case. One must be cautious in this case, though, 
because the presence of CO in the gas at all times clearly affects 
the abundance of many molecules, such as CN and N2H + . 



4 DISCUSSION AND CONCLUSIONS 

We have shown how the modelling of the depletion of heavy mole- 
cular species in dense and cold regions of molecular clouds de- 
pends on the determination of accurate adsorption energies, Ed- 
Different values of Ed affect the time molecular species reside on 
the grain surfaces and can provide a way to help explain the non- 
complete freeze-out of molecules from the gas phase. In particular, 
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1 10 2 10 4 10 6 1 10 2 10 4 10 6 

time (yr) 

Figure 1. Abundances of gas-phase (solid lines) and solid-phase (dashed lines) CO (black), NH3 (red), CN (green), CS (blue), N2H+ (purple), for the 
chemical models with T d =W K and: (a) n H = 2 X 10 4 cm" 3 ; (b) n H = 2 X 10 6 cm" 3 ; (c) n H = 2 X 10 4 cm" 3 , U2E D ; (d) n H = 2 X 10 6 cm" 3 , 
U2E D ; (e) n H = 2 X 10 4 cm" 3 , V3E D ; (/) n H = 2 X 10 6 cm" 3 , IBE D . 



we argue that the adsorption energies of molecules on a substrate 
of H2 ice can be several times lower than the ones commonly used 
in current chemical models. We also showed that H 2 could occupy 
all the available grain sites in a relatively short time, ~ 0.5 yr, 
which could produce a change in the adsorption rates of molecules 
on grain surfaces. 

We ran a few simple gas-grain models to test how much rela- 
tively small changes in the adsorption energies of molecular species 
affect the gas-phase abundances. We find that if we use values 
of Ed half of the standard values, the depletion of most of the 
molecules is greatly reduced, although it is still at levels of a fac- 
tor of 10. A further reduction of the molecular adsorption energies, 
down to a factor of 0.3-Ed, which would be a conservative value 
of how much the adsorption energies can be changed, shows that 
molecules show very little or no depletion even at relatively high 
densities, 2 x 10 6 cm -3 . Thus, lower evaporation energies of heavy 
molecular species due to the coating of grains by H2 molecules may 
explain the presence of CO and other heavy molecules in the gas in 
average dark clouds. 

As we indicated in Sect. 12.11 we do not expect to have more 



than one monolayer of H2 molecules due to the almost instanta- 
neous evaporation of any H2 landing on H2-ice (see Table [TJ. On 
the other hand, we estimate that a CO molecule landing on a grain 
completely covered by H2 will be able to scan the N s sites on the 
surface of the grain in a time, ~ 0.2 s, much shorter than the ex- 
pected evaporation time, t evap (CO) ~ 40 min. If the grain sur- 
face is only partially covered by H2-ice, the effect will be that any 
CO molecule landing on H2-ice will be able to find a free site out- 
side of the H2-ice, at which point it will remain trapped on the 
surface. This will reduce the number of available sites to be cove- 
red by H2 molecules by about 500 in a year, as we expect about 
one CO molecule landing on the surface per day at a gas density of 
10 4 cm" 3 , and should not have any effect on the basic mechanism 
we discussed. 

In the calculations shown in Sect. |2.1| we assumed a "flat" or 
homogeneous grain surface, represented by just one value of the 
absorption energy, Ed, seen by the H2 molecules. But grain sur- 
faces are thought to be rough, with an inhomogeneous structure 
of valleys, protrusions, cavities, and with and a certain degree of 
porosity. Landing H2 molecules will be then subject to a variety 
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of binding energies. Results of the effects of surface roughness on 
the formation rate o f H2 on grains fcuppen & H erbst 20051 120071: 
ICuppen et alj|2006h show that irregularities affect the binding en- 
ergies of the surface molecules, effectively increasing the binding 
energy "seen" by each molecule. We expect that H2 molecules will 
tend to fill the valleys of the surface and, probably, help to build, 
locally, more than one layer of H2 ice. 

We realise that there are several objections that could be raised 
against this approach. There were some experimental results that 
show that the su rface coverage of H2 might be ~ 20% for a binding 
energy of 450 K dGovers et alJll980l:ISchutte et aljll976h . Additio- 
nally, in the simulations of Cuppen & HerbsJ 1 2007 ), H2 molecules 
pile up with CO in it. The potential dislocation of (almost com- 
pletely) coating H2 by an incident H2 or heavier molecule can act 
against this proposed mech anism for desorption of he avy molecules 
from grains in dark clouds. iTielens & HagenNl982l) also indicated 
that the H2 molecular layer may not be rigid enough to prevent 
oncoming heavier molecules to bind themselves to the surface un- 
derneath the H2 layer. 

Nonetheless, even if the H2 layer does not completely cover 
the grain surface and cannot completely avoid that landing heavy 
molecular species bind with the surface below it, there is no doubt 
that H2 molecules must be located all over the grain surfaces, due 
to the ubiquity of H2 in molecular clouds. This should have some, 
maybe substantial, effect on the effective adsorption energy value 
of molecular species that in the right environment makes this me- 
chanism relevant. Particularly, after proving that some relatively 
small changes can affect a lot the depletion rate of heavy molecules. 

At the same time, the present mechanism and the several pre- 
viously proposed desorption mechanisms could very well work si- 
multaneously. It is yet to be seen which process dominates in a cer- 
tain environment, how effective the H2 coating can be or in which 
kind of environments it would be dominant. 

The presence of a monolayer of H2-ice on the surfaces of 
grains in molecular clouds could be observationally tested. Un- 
fortunately, there is no clear evidence yet of the dete c tion of 
solid H2 on interstellar grains. ISandford & Allamandolal dl993bl) 
showed that solid H2 could be detected by an infrared absorp- 
tion band at 2.417 fjm (4137 cm -1 ) attributed to the Qi(l) pure 
vibr ational transition of H2. This ba nd has never been detected 
and ISandford & Allamandolal dl993bl) argue that it wo uld proba- 
bly trace the H2 trapped in water ices. ISchaeferl d2007l) proposed 
six pure para-hydrogen pair transition bands as candidates to ex- 
plain unidentified emission features in ISO-SWS spectra of the 
NGC 7023 nebula. If any H2-ice absorption band was detected, it 
would be a combination of existing surface and mantle solid H2. 
Infrared spectroscopy studies of dust have found vis ual extinction 
thre sholds of sa 3 for H 2 Q ice dWhittet et alii 19881). sa 5 for CO 
ice dWhittet et al.l 19891) , and « 4 for CO? ice dWhittet et al.l2007l) . 
which roughly corre spond to very few monolayers of ice, proba- 
bly between 2 and 5 dHassel et alj2010l) . The detection of only one 
monolayer of H2 ice, in the absence of any solid H2 in the mantle, 
may not be yet feasible. 

Finally, it is interesting to consider if the H2 -coating mecha- 
nism proposed here can have any influence on the determination of 
the visual extinction thresholds discussed above. As TableQjshows, 
the evaporation time of H2 on H2-ice is highly sensitive on tempe- 
rature, t evap ~ 4.5 s at 15 K. Thus, at the expected temperatures 
of the grains at those visual extinctions thres holds (between 12 and 
20 K. lDraine & Ledll984l ; ls"mith et alj|l993l) . H 2 molecules do not 
have time to build up on the grain surfaces and we do not expect 
them to play any relevant role. Only when the grain temperatures 



go down to about 10 K will this mechanism delay the growth of 
ices 

Our simple modelling of the chemistry, changing the values of 
the adsorption energies, also proves how important is to obtain re- 
liable estimates of evaporation energies, either laboratory measure- 
ments or theoretical calculations, for heavy molecules on H2 ice or, 
for that matter, in different kinds of substrate. A simplified experi- 
ment could be the study of H2 adsorption on various types of pure 
ice (H 2 0, NH 3 , SO, SO2, H 2 CO, ...) to estimate the interaction 
potential between the H2 molecules and other heavier molecules. 
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